FLEX-description of the spectral functions near singlet-triplet transition 
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In a previous article^, we have investigated the non-equihbrium two-level Anderson model with a 
simple iterative perturbation theory. Here we use here a more sophisticated perturbative method, 
the fluctuation-exchange (FLEX) approximation. The great advantage of FLEX is its conserving 
nature, and that it can describe well the Kondo energy scale, the Kondo-temperature, Tk- As it 
was expected from the results obtained with iterative perturbation theory, the FLEX description 
can give back also the relevant features of the spectral properties. 
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I. INTRODUCTION 



In the past decades the investigation of electronic 
transport through quantum nanostructures gained a lot 
of interest both experimentally and theoretically<^"— By 
the development of nanotcchnology the properties of 
these structures gets more controllable. In most of the 
cases, the nanostructures are investigated under out of 
equilibrium conditions in the experiments. In the case of 
a quantum dot, a single atom or molecule between the 
biased leads, the system can be described well with the 
non-equilibrium Anderson modeli^^— The dynamical fea- 
tures of the quantum dot show an odd-even asymmetric 
behavior with the number of electron on the dot. If there 
are odd number of electrons residing on the dot, spin— ^ 
Kondo effect can be observed. By even occupation, on 
which we will focus in the present paper, the situation 
is more complex, the triplet Kondo efFccti^i"— and the 
Kondo-efFect associated with singlet-triplet degeneracy^^ 
were found in vertical^ and latera l-'^^'^^i^^ quantum dots 
(QDs). The lateral and vertical QDs have different ge- 
ometry, with different number of modes connecting the 
dot levels to the leads, while in the vertical case the num- 
ber of modes is more than one, in the lateral case one or 
two modes connects the QD to the conduction electrons 
of the leads. At the triplet side of singlet-triplet transi- 
tion, we can speak about fully screcne d^^i^^ (2 modes) 
and undcrscrccncd spin— 1 Kondo-cfFcct--''-'^ (1 mode), 
the transport properties are totally different in this two 
case. 

To the description of the quantum dots under out 
of equilibrium conditions, a great number of perturba- 
tive and non-pcrturbative theoretical scheme of Anderson 
model was generalized from equilibrium. The most effec- 
tive perturbative treatments are the iterative perturba- 
tive approaches which was developed also for singlei^iii 
and multi-levelii2^ Anderson models. These descriptions 
are stable even at rather high interaction parameters and 
have the great advantage that they can be easily general- 
ized to non-equilibrium using Keldysh formalism;^ but 
cannot give back precisely the width of the Kondo-peak. 
(A detailed description of the adventages and disadvent- 
ages of this method is given in RefQ.) 



There arc a lot of non-pcrturbative treatments also, 
which were generalized to out of equilibrium from suc- 
cessful equilibrium treatments. The Bethe Ansatz ap- 
proach was applied to finite biases?^ and recently it 
was worked out for the Anderson model The Nu- 
merical Rcnormalization Group technique was general- 
ized with applying a single-particle scattering basisi^ 
Other Rcnormalization Group techniques also had suc- 
cesses in the description of a steady-state non-equilibrium 
state<^ Also numerical methods such as quantum Monte 
Carlo^^ and Density Matrix RG^ was generalized to out 
of equilibrium and was applied to interacting nanostruc- 
ture transport problems, but these schemes suffer from 
the difficulties inherited from equilibrium. The power- 
ful technique of the Dynamical Mean Field theory was 
extended to non-cquilibriumi^ Also a yet another tech- 
nique, the iterative summation of real-time path integrals 
was developed^^, and its validity has been confirmed by 
a detailed comparison to approximative approaches. 

The relative simplicity of generalization to non- 
equilibrium of the perturbative schemes is also present 
in fluctuation-exchange (FLEX) approximation^ but 
FLEX is a conserving^^— approximation for the whole 
parameter space, and describes well the width of the 
Kondo peak in the equilibrium spectral function. It 
cannot describe the charge-excitations (the so-called 
Hubbard-peaks) in the spectral function^, but we as- 
sume that in the non-equilibrium properties this does not 
play an important role. FLEX sums up some classes of di- 
agrams to infinite order, we chose to sum up the particle- 
hole type diagrams, in the present case of the multi-level 
Anderson model. The other summable diagram types do 
not play a significant role in this problem. We describe 
here lateral and vertical two-level quantum dots with ap- 
proximately two electrons in average on them, while these 
systems show many features of the singlet-triplet transi- 
tion in and out of equilibrium i^^i^°i^'^i'^^i^^ 

The present paper is a second part of a two part se- 
ries, in the first pariji (called paper I in the following) 
we have investigated the same system with the scheme 
of the iterative perturbation theory (IPT) .— i^^ii^ As we 
found there, that IPT can capture many properties well 
in equilibrium, and it gives satisfactory results even in the 
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non-equilibrium limit. Based on the successes of that per- 
turbative description, we assume that also the FLEX de- 
scription can give good equilibrium and non-equilibrium 
results. 

The paper is organized as follows. In Sec. Ill Al we in- 
troduce the model of two level non-equilibrium Anderson 
model, and the vertex perturbation theory done in the 
Coulomb interaction and the Hund's rule coupling. In 
Sec. Ill BI we show how to apply the fluctuation-exchange 
approximation to this model. In Sec. Ill CI we show the 
details of the FLEX iteration loop, and in Sec. Ill Dl we 
write about the limitations of this scheme. In Sec. IIIII 
we detail the results obtained for the spectral functions, 
in Sec, nil Al for the symmetrically hybridized case, while 
in Sec. HUB] for the asymmetrically hybridized case. We 
conclude in Sec. IIVI and give some calculation details in 
Appendix and |B] 

II. THEORETICAL FRAMEWORK 
A. Model 

We describe an out of equilibrium quantum dot with 
the two-level Anderson model as we have done it in paper 
I. The Hamiltonian of the system can be splitted to a non- 
interacting and interacting term: H = Hq + Hint , where 
the non-interacting term, Hq can be further divided onto 
three terms similarly to the Eq. (1) of paper I, 

Hq = -ffcond + ^^hyb + -f^O,dot- (1) 

In this expression i/cond describes the conduction elec- 
trons of the leads, i?hyb accounts for the tunneling be- 
tween the conduction electrons and the electrons residing 
on the dot. This two part is different in the case of lateral 
and vertical quantum dot. 

For lateral QD, one electron mode is coupling to the 
dot per lead because the lateral QD is close to pinch-off, 
the Hamiltonians are (cf. Eqs. (3)- (4) of paper I), 

-^cond = ^ ] ?Qc|Q,£rC{Q(T 1 (2) 

Hl:^b - E ^c.^{cl^A- + h.c.) . (3) 

c|^g. is the creation operator of a conduction electron 
in the left or right lead {a G {L,R)) with energy ^q, ~ 
C + Mai i^^a = eVa IS the bias applied on lead a, the non- 
equilibrium condition is taken into account in a different 
shift of the chemical potential in each lead) and spin a. 
dia destroys an electron with spin a on the dot level i G 
(-f , — ) and tai is the tunneling matrix element connecting 
the lead a and dot level i, i is assumed to be spin- and 
energy independent. 

In the vertical case, the dot is connected to the leads 
with a relatively large surface, therefore there are many 
different conduction channels coupling to a given dot 



state. For each dot state i and lead a one can construct 
a simple linear combination of the modes, c| j^^, that 
hybridizes with dia, and assumed to be independent for 
i = ±. The Hamiltonians are (cf. Eq. (5)-(6) of paper 

I), 

^cond — £,cxC^^io,fjC^ioi(T , (4) 

^hyb = E tc.^{clo.A^ + h.C.) . (5) 

Eqs. ([2]) and (|4]) shows that the interaction between the 
conduction electrons is neglected. 

The last term of Eq. ([T|) describe the bare dot levels, 

^^o.dot = E^y'^L'^j'^ . (6) 

The eigenvalues of the e are the energies corresponding 
to i = ±. 

Throughout this paper we consider a special com- 
pletely symmetric tunneling matrix, similarly to the IPT 
description of paper I, because this case captures most of 
the spectral and transport properties obtained in exper- 
iments for a two-level QD. We assume, that one of the 
dot levels is even the other is odd (— ) under reflec- 
tion. Under these conditions, the tunneling matrix has 
simple structure: t^^^ = tii,+ and — — In 
one of the cases we specialize further this situation with 
\ta,+ \ = \ta,- \ which gives a particle-hole symmetric case 
with exactly two electrons on the dot. By these values of 
the parameters, we define here the hybridization param- 
eter, the F-matrix, connected to the level width of the 
dot states, 

r., = 2^p(0)Ei-C, . (7) 

a 

where p(0) is the spectral density of the conduction elec- 
trons at Fermi-energy, we take it p[Q) = 1. Here is given 
a general form, but we have to note that by the above 
detailed special parameter values, only the diagonal el- 
ements are finite of the matrix. The diagonal elements 
will be called, for the sake of brevity, r+ = r++ and 

r„ = r__. 

Using the expression (7) of paper I, the electron- 
electron interaction can described with 




where the occupation number is riio- = dj^dia and the 
spin is = i ^ dl^a^^idi^i with the Pauli matrices, a ~ 

iacr' 

(ax,cFy,<Jz)- U is the on-site Coulomb- interaction, while 
J > denotes the ferromagnetic Hund's rule coupling. 
Eq. (jH) contain explicitly that the relevant case of this 
description is that, when there is (approximately) two 
electrons on the dot. 




FIG. 1: The construction of the particle-hole like vertex used 
by building FLEX self-energies 



To introduce a systematic perturbation theory, which 
treats the on-site Coulomb interaction and Hund's rule 
coupling on equal footing, we write Eq. to normal 
ordered form, 
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(9) 



where the normal order part will play the role of a per- 
turbation, the second term is a simple renormalization of 
the single electron levels. 

The normal ordered term of Eq. Q can be rewritten 
using an antisymmetrized interaction vertex of particle- 



hole type, F;' 



shown in Fig. [TJ as 
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which is very similar to the expression (11) of paper I, 
but in opposition to that, this expression uses particle- 
hole vertices, instead of particle-particle vertices, i'l^ 
is antisymmetrical under exchanges ia ^ ma' and ja' -O- 
na, the finite matrix elements of the interaction vertex 
can be easily found. Creating the diagrams with this 
vertex function, the details of the on-site Coulomb in- 
teraction and the Hund's rule coupling will be hidden in 
the vertex (it can be resolved with the substitution of 
the diagrammatic form of the vertex in Fig. [T]), and the 
equal-footing treatment of U and J is ensured. 



FIG. 2: Diagrammatic representation of Il-function 



are relevant, therefore, we consider only this group of dia- 
grams in the following. We can build equivalent diagrams 
with the two kind of interaction vertices which was used 
in paper I and that in the recent work in Fig. [T] In op- 
position with the particle-particle type vertex of paper I, 
the particle-hole type behaves f = — f by exchange 
of the two outgoing legs. Where a, /?, 7, 6 are compound 
indices, a = (/, a, k). The vertex has simple behavior in 
Keldysh-space: 
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(11) 



While the interaction vertices are spin-conserving, only 
vertex elements with special spin arrangement are finite. 
We arrange the vertex elements in a matrix based on the 
spin arrangements. The rows (upper pair of spins in the 
vertex) and the columns (lower pair of spins in the vertex) 
have the spins ttj ti: it and 4,4. respectively. We empha- 
size that in each position of this spin matrix, a matrix 
of the level and Keldysh indices resides. Because of the 
special spin arrangements this matrix can be described 
by two independent matrix variables as follows. 
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(12) 



If we substitute in the finite vertex elements, we can find 
the finite elements of these singlet and triplet matrices 
occurring in the matrix 



(13) 
(14) 
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B. Out of equilibrium fluctuation exchange 
approximation 

In this section we will show how to apply the 
fluctuation-exchange approximation (FLEX) to the de- 
scription of a two-level QD under non-equilibrium condi- 
tions. It is shown how to calculate the out of equilibrium 
Green's functions {G'^''^ , k and k' sign the branches of 
the Keldysh-contour) with FLEX-approximation in the 
self energy. 

The great advantage of FLEX that it sums up a group 
of diagrams to the infinite order. In the case of the two- 
level Anderson model the so-called particle-hole diagrams 



With this classification of the spin, the whole treatment 
simplifies. Calculating with this singlet-triplet matrices 
we have to treat only Keldysh and level indices in the 
calculations. 

By the construction of the diagrams, it is convenient 
to introduce a new quantity, the so-called polarization 
function, defined diagrammatically in Fig. [2l 



^(°)::S(^-t') = ^^G-(^-^')G^^(t'-^) 



(15) 



where G^J is the full Green's function. Using the po- 
larization function and the interaction vertex, it is very 
straightforward to construct the full FLEX self-energy, 
shown in Fig. [31 In the figure the hierarchical order of 
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FIG. 3: FLEX self-energy 
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FIG. 4: Full self-energy 



the self-energy can be observed. It is relevant to note 
that using these definitions, the second order diagram is 
over-counted in the FLEX self-energy, this will have to 
be compensated in the full self-energy, as shown in the 
diagrammatic form of it in Fig. |3] (all the other orders of 
the FLEX self-energy are not over-counted). 

The n -\- 1-th order diagram of the FLEX self-energy 
can be described in the time-domain 

oo 



X n(")t'k_,(t„-t')fff,„Gi(t-t'), (16) 

where the integration goes over the inner variables, 
ti . . .tn- The part of this expression which depends on 
the inner variables can be organized in a new quantity 
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(17) 



diagrammatically shown in Fig. [5l Using this expression, 
we can write down n + 1th order self-energy in a simpler 
form, 

'^a.FLEX^^ ^ ~ 

= -Cf°n(")^;|;(t-Of£^„Gf(i-i'). (18) 



The polarization functions have the same structure in 
spin as the interaction vertices, therefore we can use the 
same matrix description for them. 
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Considering only the different inner part described by 
in each order, it is transparent that it can be 
summed up easily to the infinite order, because they give 
a geometrical order: 



n(w) 



n(") = n(°)(w) + n(°)(w)(-if)n(°)(w) 
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n(°)(L.) fi + ifn(")(o 



(20) 



Using the singlet-triplet description in the polarization 
matrix (H^, n^), similarly to the interaction vertices, 
Eq. (PO)) reduces to the two equations as follows. 
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(21) 
(22) 



The inner structure of 11 and F matrices can be 

=s , t =s , t 

described by a compound index, which contain level and 
Keldysh indices. If we describe 11;^ as a 2-index matrix, 
the compound indices are 
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where I = ± are the level indices, k = 1, 2 arc the Keldysh 
indices. 

While the inner part brings the order of the diagram, 
we can obtain the full FLEX self-energy by summing up 
these expressions to infinite order. 



^FLEX 



(n+l) 
FLEX ' 



(23) 



n=l 



V/3 

^qFLEX 



(t-f) = -KrKi'k(t-t')x (24) 



(25) 



Details of the FLEX iteration 



FIG. 5: nth order polarization function, II'"' 



In a numerical calculation the Green's functions are 
represented on a finite, discrete frequency mesh charac- 
terized by two basic quantities, the number of the grid 
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points (N) and the width of the frequency range the func- 
tion was represented on (Q). We used an equidistant 
mesh with distance between the grid points Aui = fl/N. 
In the present calculations we use f2 = 1000 and N = 2^^ 
to obtain spectral functions. 

Here we write down in detail one iteration loop of the 
FLEX approximation. It is relevant to emphasize in the 
beginning of this section, that the full Green's function is 
iterated in the FLEX scheme, every points of the function 
have to converge. We use the w-space Green's functions 
to check the convergence. As a convergence criterion, we 
are searching for the largest difference of the previous and 
the new Green's functions as a function of w-points. The 
initialization of the full Green's function at the beginning 
of each step 



where the components are 



(26) 
(27) 



where the number of the iteration signed in the upper 
index. In the 0th order we initialize the Green's function 
with the hybridized Green's function, the calculation of 
it is introduced in the Appendix El To define the po- 
larization we need Green's functions in 'i'-space, using 
Fourier transformation: 



(28) 



Using a flip, we calculate G"^J(i) ^ G"°l{-t), and the 
polarization in time-space is obtained from Eq. (jl5|) . In 
the following, we will use matrix notation for the po- 
larization and we will need a singlet and triplet part of 
n^°\ which will have only level and Keldysh indices. The 
Green's functions are diagonal in spin, therefore by def- 
inition, the 0th order polarization has the same value in 
the singlet and triplet channel as follows. 



2(0) 
n(o) 



(29) 
(30) 
(31) 



To obtain the interacting polarization function we need 
R^liuj) = FFT{^fl{t)), and using the above defined £^ 

and r^, we can create the fully interacting n following 
Eqs.lni) and dm). 

' construct the components of the self-energy the po- 
larization function is needed in time-space. 



(32) 



using this function, the full self-energy can be constructed 
from four parts 



^ — ' \ — / aa 



(34) 



a/3 



^un-m = -[^ + ^) Si^)5^,Sit) , (35) 



^'Kit) = S..Y.C^l (36) 

al3 

^'<(t) = -5:(rn(°)(t)r)''Gf(0, (37) 



where Eq. ([M)) describes the FLEX self-energy summed 
up to infinite order, Eq. (j35p contains the linear shift 
obtained by taking the interaction to normal order, Eq. 
(1551) describes the first order Hartrec diagrams, and Eq. 
(1371) the second order diagram, which compensate the 
over-counting. The FLEX ((Ml), the first order (|36l) and 
the second order ([57)1 self-energies can be calculated in a 
simple way using the singlet-triplet representation of the 
interaction vertex and polarization matrix. This proce- 
dure is shown in the Appendix [BJ 

The Dyson equation needs the total self-energy in the 
'w'-space. 



J:''iiuj) = FFT (S"^(t)) 



(38) 



as follows. 



G-'"^'!(^)=5-'a(^)-S"^(^), (39) 

where we used the hybridized Green's function (g). With 

an inversion, we can have G"+if(w). Calculating the 
convergence parameter in the nth loop, 

^\G-+^iiu^)-G-iiu^W^ . (40) 

At the end, we interpolate to have the Green's function 
for the next iteration. 



aG"+if(c.) + (l-a)G"^(a.), (41) 



In the iteration procedure we iterate this loop until the 
e" become small enough, in some nth order, which means 
that the Green's function is converged. When we reach 
this point we can extract spectral functions from the re- 
sulting Green's function from 



-ilm GlM, 



(42) 



S"^(t) = J^FLExZW + ^ct'-'ait) 
+ S<(t)-S<(t), 



^1 



(33) 



in the followings in many of the cases, we will use the to- 
tal spectral function, /Ot(w) = a- Pio-{^)i or a spectral 
function only summed up for the spins. 
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D. Limitations 

We have seen in paper I that the IPT suffers from insta- 
bihty at large values of U /T and J/F parameters. Similar 
kind of instability exists in the FLEX scheme: from ar- 
bitrary starting Green's function the iteration converges 
only at small values of U /F and J/F. The error caused 
by the not satisfactory behavior of the self-energy con- 
tributions for large U /T and J/F parameter rates, if we 
use the hybridized Green's function to obtain the self- 
energy. But in opposition to IPT, starting from a con- 
verged Green's function, high enough interaction param- 
eters can be reached by increasing the interaction pa- 
rameter in small steps. This means that approximately 
U/T K, 2 can be reach using this procedure starting from 
a Green's function converged at U/T k, 0.75. Also the 
J/F rate has to be small enough to have a good conver- 
gence for the starting Green's function. 

The unit of energy will be the on-site Coulomb inter- 
action U . For the present results of this paper, the cal- 
culations have the following two part. In the first part, 
we calculate a FLEX Green's functions starting from the 
hybridization Green's function corresponding to the pa- 
rameter set with U = 1 and with comparatively high 
T± values to reach conductance from this somehow arbi- 
trary Green's function which is far from the final FLEX 
Green's function. After we have the converged FLEX 
Green's function we can improve the U IT± ratio with 
increasing the value of U . We increased U step by step 
with 5U = 0.1 and the final value was U = 2. 




FIG. 6: The three singlet and three triplet states which can 
be constructed on the dot. The singlet states are labeled with 
the occupied levels, the triplet states with the Sz value of the 
spin. 



sinjilcl 



-> 







A-2J 



FIG. 7: Different regimes corresponding to different ground 
state configurations as a function of A for fixed J 



III. RESULTS AND DISCUSSION 

When the levels are chosen to be particle-hole sym- 
metric, i. e. e± = —1.5U ± A/2, = e ^ = 

and \ta+\ = \ta-\, two electrons are residing on the dot, 
or when \ta+\ ^ \ta-\, the number of electrons is close 
to two, and the states can be described well with ef- 
fectively two electrons. These two electrons can form 
three 5* = singlet-state and three 5* = 1 triplet state 
on the dot. The singlet states can be represented in 
terms of electron creation operators as |s_) = d^_^d^_^\0), 

\s+) = <xjo), = ^K^dL^ - di^di^m, 

indicating the involved levels by the subscripts, while 
the triplet states can be written as |t, 1) = d^^d^^jO), 

ii,o) ^ -i=(4/_^ + jt_^4^)io), |t,-i) = 4/_,io), 

labeled by the 5*2 quantum number and the vacuum, |0), 
refers to the state where the investigated two levels are 
not occupied. 

The energies of these states can be distinguished by 
the values of the Hund's rule coupling (J) and the level- 
splitting (A), while the on-site Coulomb U gives equal 
contribution to all of the energies of these states. In 
Fig. iniwe can see the evolution of the energy-levels as a 
function of the Hund's rule coupling. Only the energy 
of the triplet states depend on J, because the energy 



term connected to Hund's rule coupling is —JS{S + 1), 
so it decreases the energy in the case of parallel spin 
alignment, S = 1, and has no contribution when S = 0. 
By the investigation of the spectral functions, we will 
change the value of A in the figures, while in the exper- 
iments the level splitting can be shifted by varying the 
parametcrsi^"— It can be seen in Fig. |6]that A/2 = J 
is a special point considering the ground state. When 
J < A/2, the ground state of the dot is non-degenerate 
singlet, but for J > A/2 the threefold degenerate triplet 
gives the ground state. For J = A/2 the |s_) and the 
three triplet states {\t,m)) are degenerate and therefore 
give a fourfold degenerate ground state. This is exactly 
true only for a bare dot, by the introduction of the hy- 
bridization, the transition between the ground states not 
occurs at the exact point of A/2 = J, but takes place in a 
wider range of the A parameter and becomes a crossover. 
The introduction of the hybridization leads the system to 
show Kondo-effect in the case of degenerate ground state. 

Fig.[7]shows what regimes can be observed in the quan- 
tum dot coupled to two conduction modes for different 
parameters. In this figure, the relevant parameter is 
A — 2J, while we take J fixed and change A. If we 
compare it with Fig. [5] we can see that for a bare dot in 
the case of < A — 2 J the ground state is a singlet, and 
the threefold degenerate triplet for > A — 2 J. Fig. [7] 
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shows an intermediate regime also, where |s_) and the 
three triplet states are degenerate, this would be the case 
only in the point 2 J = A for the bare dot, but for finite 
hybridization, this regime becomes finite. 

The first regime of Fig.[7]with triplet ground state cor- 
responds to triplet Kondo effect. In this case there is a 
finite S = 1 spin on the dot, but it is screened by the two 
conduction electron modes connected to the dot. Each 
mode screens S — ^ from the whole spin, i. e. two sepa- 
rate Kondo effect screens the whole spin of the dot, which 
can be characterized by different Kondo-temperatures, 
and because the hybridization-parameters, T±, 
defined through Eq. ([7]) can be different for the two levels. 
This phenomena is called two-stage Kondo effect J^i^Si^ 

The second regime, connected to the singlet-triplet 
Kondo effect in Fig. [7] corresponds to the case where 
|s_) and the three triplet states, \t,m) m = —1,0,1, 
gives an effectively degenerate ground state. They can 
be degenerate in a finite range of A because of the fi- 
nite value of T±. The Kondo-temperature correspond- 
ing to this situation, T^, gives the energy scale of the 
singlet-triplet Kondo-effect. Generally is larger than 
T^i^ The system shows Kondo-effect for the parameters, 
1 2 J — A I < T^, which means considering the effect of the 
hybridization, that 2J is effectively equal to A, the four 
above mentioned states can be considered degenerate. 

The third regime, called singlet, is corresponding to the 
case, where |s_) is the ground state, when |2J — A| > 
it cannot be effectively degenerate with any other states. 
The first excited state, which is the threefold degenerate 
triplet, can be reached by applying finite energy, this has 
a clear fingerprint in the spectral functions. 

A. The symmetric case (r+ — r_) 

In this section we assume that the dot is completely 
symmetric, the diagonal elements of the F-matrix defined 
in Eq. ([7]) are equal, which means that the tunneling 
matrix has the following structure. 



We investigate the spectral functions for J/U = and 
J/U = 0.15. In both of the cases the hybridization 
parameters are F±/J7 = 0.785. It is a clear feature of 
the FLEX spectral functions, that the FLEX approxi- 
mation cannot give back the charge excitation peaks^^^ 
similarly to all of the approximations which calculate 
the end Green's function by iterating the fully interact- 
ing Green's function in every stepi^"— Otherwise such 
treatments are the only candidates to be conserving, we 
assume that the above mentioned feature doesn't influ- 
ence the description of the transport properties too much. 
Beside this feature, the FLEX approximation describes 
well the width of the Kondo-peak, which is connected 
to the Kondo energy scale and what is assumed to be 
relevant by the description of the transport properties. 




FIG. 8: Total spectral function, pri^^) for J/U = and 
r±/U = 0.785 for different values of level splitting, A/U 

Considering first the spectral functions of Fig. [S] for 
J/U = and different A/U parameters. We have to 
note at the beginning that in the lack of finite J, only 
two regimes can be investigated from the three regimes 
of Fig. [3 as we change A, the triplet Kondo effect is 
not present for J/U = 0. For A/U = and small level 
splitting values (A < T^), the spectral functions show 
the existence of a singlet-triplet Kondo effect, while they 
show a one-peak structure with a peak at the Fermi- 
energy, which is the Kondo-peak. Above some value of 
A/U, the Kondo-peak splits. For large A/U values the 
two side-peaks correspond to the triplet excitations, and 
they are residing on the exact place of w = ±A/2. With 
this finite amount of energy we can excite the first ex- 
cited state which is the threefold degenerate triplet in 
this case. For intermediate A/U values (A/J7 = 0.6 in 
Fig. [8|) the splitting is not corresponding unambiguously 
to the triplet excitations, the peak positions are influ- 
enced by the Kondo effect also, while the system is in 
the transition regime between the singlet-triplet Kondo 
and singlet regimes. 

In the presence of finite Hund's rule coupling, J. the 
spectral function is shown on Fig. [9] for J/U ~ 0.15. 
In this set of spectral functions we can see differences 
compared to Fig. [51 By the introduction a finite value of 
J, we can investigate the triplet Kondo-effect for << 
2J- A (r+ = for r+ = F_), which is true for smaU 
values of A for the spectral functions of Fig. [9l Which 
was assumed earlier, that is smaller than T^, can 
be clearly seen for A/[/ = in Fig. [HI i.e. the width of 
the triplet Kondo-peak is smaller than the width of the 
singlet-triplet Kondo peak (the effect is larger for larger 
values of J). For large values of A (e.g. A/U = 1.2 
in Fig. in]), we can see the above detailed fingerprints of 
singlet ground state, the two triplet side peaks at uj/U = 
±(A — 2J)/U. For intermediate A values the spectral 
function is influenced by the singlet-triplet Kondo-effect, 
when \2J - A| < T^. For A/U = 0.6 this is the case, 
and the corresponding spectral function on Fig. ^ has a 
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wider Kondo-peak than for A/[/ = 0, as it was assuraed. nL i I i I i I i 

-2 -10 12 

co/[/ 



B. The asymmetric case, r+ 7^ r_ 

In the case r_|- 7^ r_, there is no particle- hole symme- 
try in the system. The vahie of the total occupation of the 
quantum dot is not exactly two any more, and the spec- 
tral function is not symmetric. For r+ ^ r_ we consider 
only the physically interesting case, when the Hund's rule 
coupling is finite. The spectral functions separately for 
the levels for parameters r+/C/ = 1.57, r_/C/ = 0.785 
and J jU = 0.15 can be seen in Fig. [TUk while in Fig.llOb 
the total spectral functions is observable. For A/t/ = 
similarly to the r+ — r_ case of Fig. [HI we can see the 
fingerprints of a triplet Kondo effect. If we consider the 
A/[/ = spectral functions for each level in Fig. [TOk. we 
can see the two different width of the Kondo-peaks due 
to the different Kondo-temperatures. This feature shows 
the dependence of on the hybridization parameters, 
F-i-. For small splitting, this is the cause of the asymme- 
try of the total spectral function, but while the triplet ex- 
citation peaks also have different width, the asymmetry 
is characteristic for all of A/C/ values. The trends with 
increasing A can be clearly observed on the separately 
shown spectral functions of Fig. [TUk . It is transparent 
in that figure, that when the Kondo-peak changes little 
and it leaves the Fermi-energy slightly then the system 
shows triplet or singlet-triplet Kondo-effect. At the A 
value where the Kondo-peak gets wider, that is the point 
of transition from triplet to singlet-triplet Kondo effect. 
For A/[/ = 1.2 where the peak of the spectral functions 
in Fig. [TUb far from w = and from the slightly splitted 
peaks, and have larger width than the peaks correspond- 
ing to Kondo effect, we can assume, that these peaks 
corresponds to the triplet excitations, and the ground 
state of the quantum dot is a singlet for this parameter 
set. 



FIG. 10: The spectral functions for the levels separately 
(a), and the total spectral function, /9t(i^) for J jU = 0.15, 
r+/[/ = 1.57 and V-jU = 0.785 for different values of level 
splitting, A/f/ 



IV. CONCLUSIONS 

In this paper we applied fluctuation-exchange approx- 
imation to the non-equilibrium Anderson model. We 
generalized the FLEX approximation using Keldysh the- 
ory. The spectral functions were investigated for different 
parameter sets, we could catch singlet-triplet transition 
by varying A parameter, and could observe triplet and 
singlet-triplet Kondo effect on the spectral functions. It 
is relevant to note as an evidence of the good description 
of the Kondo energy scales, that the width of the spec- 
tral functions corresponding to triplet and singlet-triplet 
Kondo effect are distinguishable. The good description of 
the width of the peaks is also relevant by distinguishing 
the triplet excitation peaks from the peaks correspond- 
ing to singlet-triplet Kondo effect. We can conclude that 
although the FLEX approximation does not give back 
the so-called charge side-peaks in the spectral functions, 
the description of the behavior the central Kondo-peak 
is satisfactory for all parameters. 
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Appendix A: Calculation of the hybridized Green's 
function 



We can construct the hybridized Green's functions in 
some easy steps. The inverse of the atomic Green's func- 
tion is (independent from spin) 



JO) 



aK (^) = i^^^J - £y) ^kk'S'(k) , (Al) 



where S{k) is the Keldysh sign. The hybridization self- 
energy is (also independent from spin) 



aeL.R 



(A2) 



where 



S++(c.) = m{2Uuj) - 1) , (A3) 

S]+-(c^) = -2mUuj), (A4) 

I]-+(c^) = -2^z(/„(w) - 1) , (A5) 

I]--(c.) = 7r*(2/,(w) ™ 1) , (A6) 

where /^(a;) = f{uj — fia) is the Fermi function. The 
hybridized Green's function is 



g 



)) '5kk'»S'(k) ^ ^ taitaj^a 



aeL,R 



independent and therefore diagonal in 'cr'. The indices 
'jctk' will be taken to a compound index 'a'. 



Appendix B: Simplification of the self-energy 
calculation 



Here we show how can we use the simple structure of 
the polarization and interaction vertex matrix to some- 
how simplify the self-energy expressions. We can use a 



fact that multiplication of 11 and F will also have 

^ =s,t =s,t 

singlet or triplet structure. The simplifications shown 
here can be used by each term of self energy which con- 
tain polarization and/or vertex matrices, we will show 
the simplification step in the FLEX self energy, 

^FLEX^it) = -J2(tUmY' Gi(t) . (Bl) 
— \ / eta 

a/3 

Because the self-energy is independent of the spin, we 
are free to fix spin e. g. in Only considering the spin 
arrangements in the above expression 



A = - (£fl(^)£)|| + {mm)l\ Gi , (B2) 



the multiplications in the brackets can be expressed by 
the singlet and triplet terms of the multiplication, 

(£S(i)£) = (£S(i)£) ■ (B4) 

The Green's functions are not spin dependent so we can 
simplify the original expression of the self energy (jBl[) 
by a form which uses the easily computable singlet and 
triplet contributions of the polarization and vertex ma- 
trices. 



S_FL£;Xa(0 



/3/3 



a/3 



This means that we have to work only with the level 
and Keldysh indices by the calculation of the self-energy, 
which acts as a very large simplification by the storage 
of the matrices. 
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